Pseudokarst topography in a humid environment caused by contaminant-induced colloidal dispersion

Over fifty small sinkholes (~1 meter in depth and width) were found in conjunction with structural damage to homes in an area south of Cleveland, TX. The local geology lacks carbonate and evaporite deposits associated with normal sinkhole development through dissolution. The morphology and distribution of sinkholes, and the geologic setting of the site are consistent with piping erosion. However, the site lacked the significant hydraulic gradient or exit points for sediment associated with traditional piping erosion. In areas of sinkholes, geophysical measurements of apparent electrical conductivity delineated anomalously high conductivity levels that are interpreted as a brine release from a nearby oil-field waste injection well. The contaminated areas have sodium adsorption ratios (SAR) as high as 19, compared to background levels of 3. Sodium has been shown to cause dispersion of soil colloids, allowing for sediment transport at very low velocities. Thus, subsurface erosion of dispersed sediment could be possible without significant hydraulic gradients. This hypothesis is backed by the observation of the depletion of colloidal particles within the E-horizon of sinkholes. However, there is a lack of precedence of waste brines initiating colloid dispersion. Also, sodium dispersion is not thought to be an important process in piping erosion in humid settings such as this one. Therefore, laboratory experiments on samples from the site area, designed to simulate field conditions, were conducted to measure dispersion verses pH, SAR and electrical conductivity (EC). Analysis of the experimental data with neural networks showed that an increase in SAR did increase dispersion. A dispersion prediction map, constructed with the trained neural network and calibrated geophysical data, showed correlation between sinkhole locations and increased predicted dispersion. This research indicates that a contaminant high in sodium content has caused colloidal dispersion, which may have allowed nontraditional subsurface erosion to occur in an area lacking a significant hydraulic gradient

GPR Method for the Detection and Characterization of Fractures and Karst Features: Polarimetry, Attribute Extraction, Inverse Modeling and Data Mining Techniques

The presence of fractures, joints and karst features within rock strongly influence the hydraulic and mechanical behavior of a rock mass, and there is a strong desire to characterize these features in a noninvasive manner, such as by using ground penetrating radar (GPR). These features can alter the incident waveform and polarization of the GPR signal depending on the aperture, fill and orientation of the features. The GPR methods developed here focus on changes in waveform, polarization or texture that can improve the detection and discrimination of these features within rock bodies. These new methods are utilized to better understand the interaction of an invasive shrub, Juniperus ashei, with subsurface flow conduits at an ecohydrologic experimentation plot situated on the limestone of the Edwards Aquifer, central Texas. First, a coherency algorithm is developed for polarimetric GPR that uses the largest eigenvalue of a scattering matrix in the calculation of coherence. This coherency is sensitive to waveshape and unbiased by the polarization of the GPR antennas, and it shows improvement over scalar coherency in detection of possible conduits in the plot data. Second, a method is described for full-waveform inversion of transmission data to quantitatively determine fracture aperture and electromagnetic properties of the fill, based on a thin-layer model. This inversion method is validated on synthetic data, and the results from field data at the experimentation plot show consistency with the reflection data. Finally, growing hierarchical self-organizing maps (GHSOM) are applied to the GPR data to discover new patterns indicative of subsurface features, without representative examples. The GHSOMs are able to distinguish patterns indicating soil filled cavities within the limestone. Using these methods, locations of soil filled cavities and the dominant flow conduits were indentified. This information helps to reconcile previous hydrologic experiments conducted at the site. Additionally, the GPR and hydrologic experiments suggests that Juniperus ashei significantly impacts infiltration by redirecting flow towards its roots occupying conduits and soil bodies within the rock. This research demonstrates that GPR provides a noninvasive tool that can improve future subsurface experimentation